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Abstract 

A general model independent approach using the 'off-shell Bethe An- 
satz' is presented to obtain an integral representation of generalized form 
factors. The general techniques are applied to the quantum sine-Gordon 
model alias the massive Thirring model. Exact expressions of all matrix 
elements are obtained for several local operators. In particular soliton 
form factors of charge-less operators as for example all higher currents are 
investigated. It turns out that the various local operators correspond to 
specific scalar functions called p-functions. The identification of the local 
operators is performed. In particular the exact results are checked with 
Feynman graph expansion and full agreement is found. Furthermore all 
eigenvalues of the infinitely many conserved charges are calculated and 
the results agree with what is expected from the classical case. Within 
the frame work of integrable quantum field theories a general model in- 
dependent 'crossing' formula is derived. Furthermore the 'bound state 
intertwiners' are introduced and the bound state form factors are inves- 
tigated. The general results are again applied to the sine-Gordon model. 
The integrations are performed and in particular for the lowest breathers 
a simple formula for generalized form factors is obtained. 

PACS: ll.10.-z; ll.10.Kk; 11.55.Ds 

Keywords: Integrable quantum field theory, Form factors 



'Permanent address: Yerevan Physics Institute, Alikhanian Brothers 2, Yerevan, 375036 
Armenia. 

^e-mail: babujian@lx2.yerphi.am, babujian@physik.fu-berlin.de 
t e-mail: karowski@physik.fu-berlin.de 



1 



1 Introduction 



Form factors for integrable model in 1+1 dimensions were first investigated by 
Vergeles and Gryanik p] for the sinh-Gordon model and by Weisz for the 
sine- Gordon model. The 'form factor program' was formulated in || where the 
concept of generalized form factors was introduced. In that article consistency 
equations were formulated which are expected to be satisfied by these objects. 
Thereafter this approach was developed further and studied in the context of 
several explicit models by Smirnov who proposed the form factor equations 
(i) — (v) (see below) as extensions of similar formulae in the original article ||. 
The formulae were proven in [^J. In the last decade a large number of articles 
were published on form factors (see e.g. references in [^)). More recent papers 
on solitonic matrix elements in the sine-Gordon model are || ^1 and for the 
<Si7(2)-Thirring model || ffl, [HJ. Also there is a nice application [H], |l2| of form 
factors in condensed matter physics. The one dimensional Mott insulators can 
be described in terms of the quantum sine-Gordon model. 

In the present article the new approach to the 'form factor program' pre- 
sented in Q is developed further. It uses the 'off-shell Bethe Ansatz' to obtain 
an integral representation of generalized form factors. The approach applies 
also to general integrable models in 1+1 dimensions where the nested version 
[ p"3[ of the 'off-shell Bethe Ansatz has to be used. Applications of the general 
case will be published elsewhere [jli], [l5|, here we restrict ourselves essentially 
to the simple no-nested version. That means the general techniques are ap- 
plied to the quantum sine-Gordon alias massive Thirring model. The article 
is a continuation of the previous one Q where the soliton field (an operator 
with nonvanishing charge) has been investigated. Here exact expressions of all 
matrix elements are obtained for several charge- less local operators. 

We repeat the investigation of the current and the energy momentum tensor 
which have been discussed before by SmirnovQ].[] Our main new results are: 
1) We propose the form factors of the local operators ipip(x), ipj 5 ip(x) and of 
the infinitely many conserved currents. 2) In order to identify the operators we 
perform Fcynman graph expansions and compare the results with expansions 
of the exact expressions. 3) We calculate all eigenvalues of the infinitely many 
conserved charges. 4) Within the frame work of integrable quantum field theo- 
ries we derive a general model independent 'crossing' formula (correcting for a 
sign mistake in a similar formula proposed by Smirnov[Q). 5) We develop the 
concept of the 'bound state intertwiner' in the context of bound state form fac- 
tors and prove several relations. 6) Using these techniques we derive a formula 
for breather form factors. 

The 'form factor program' is part of 'bootstrap program' for integrable quan- 
tum field theories in 1+1-dimensions. This program classifies integrable quan- 
tum field theoretic models and in addition it provides their explicit exact so- 

1 In this work similar integral representations of form factors were proposed. We have 
checked that for small coupling the results for the four-particle matrix element of the current 
and the energy momentum tensor agree with ours. We could not prove that both representa- 
tions agree in general. 
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lutions in term of all Wightman functions. These results are obtained in three 
steps: 



1. The S-matrix is calculated by means of general properties as unitarity 
and crossing, the Yang-Baxter equations (which are a consequence of in- 
tegrability) and the additional assumption of 'maximal analyticity'. This 
means that the two-particle S-matrix is an analytic function in the physi- 
cal plane (of the Mandelstam variable (P1+P2) 2 ) and possesses there only 
those poles which are of physical origin. 

2. Generalized form factors which are matrix elements of local operators 

out ( P ' m ,...,p' 1 \0(x)\ Pl ,...,p n ) in 

are calculated using the S-matrix obtained in 1. More precisely, the equa- 
tions (i) — (v) given below on page || are used as solved. These equa- 
tions follow from LSZ-assumptions and again the additional assumption 
of 'maximal analyticity' (see also ||). 

3. The Wightman functions are obtained by inserting a complete set of in- 
termediate states. In particular the two point function for an hermitian 
operator 0(x) reads 



;o| OWO (o,,o,^I/.../ 7 ,|L_ r * k _ 



x (0\G(0)\pi,..., Pn ) m e-^P*. (1) 

The on-shell program i.e. the exact determination of the scattering matrix 
using the Yang-Baxter equation was formulated in Jig, |l7 | (for reviews see also 
[ |l8[ [Tg| ] ) . Off-shell considerations were carried out in |J, |2j] and in j| [2(J , where 
the concept of a generalized form factor was introduced. The explicit evaluation 
of all the integrals and the sum in (|l|) remains an open challenge. A progress 
towards a solution of this problem has recently been achieved by Korepin et al. 
[ pl| . Up to now it has even not been proven that the sum over the intermediate 
states converges.^ We expect that our new representations of form factors will 
help to solve these problems. The 'bootstrap program' does not use classical 
Lagrangians and any quantization procedure to construct the quantum models. 
We have contact with the classical models only, when at the end we compare 
our exact results with Feynman graph expansions which are based on these 
Lagrangians. 

In the previous paper [^) an integral representation for general soliton matrix 
elements of the fundamental fermi-field of the massive Thirring model has been 
proposed. In the present paper we generalize this formula and investigate in 
particular charge-less local operators. The strategy is as follows: 

2 However, it is known that the higher particle contributions are very small compared 
to the leading ones. 
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For a state of n particles of kind oti with momenta pi — m sinh and a local 
operator 0(x) the generalized form factor is defined by 



;0|O(a0|ai(pi), 



a n ( Pn ) ) m = e - fa ^+-+P»)O a (g) , for 9 1 > ■ ■ ■ > 6 n . 

9i, . . . , n ) has been used. 



where the short notation a = (a\, . . . , a n ) and 
We make the Ansatz 



O«(0 



Co 



dz m h{tz)p°{e_,z) *„(6,z) 



Co 



with the Bethe state ^aidi?) defined by eq. ( p^ ) and the integration contours 
Cg_ of figure [l]. The scalar function h(9,z) (see eqs. (p^)-(|l7|)) is uniquely deter- 
mined by the S-matrix and the 'p- function' p°(0_,z_) depends on the operator 
0(x). By means of the Ansatz we transform the properties (i) — (v) of the 
co-vector valued function Oai&) to properties (i 1 ) — (V) of the scalar function 
p°(6_,z) which are easily solved. In particular we obtain the p-functions for the 
local operators^ Af ffiip] (x), Af \jjj"f 5 ip] (x), the current j M (x) = J\f [^7 M V] ( x ) 



the energy momentum tensor T^ u {x) — £Af 



infinitely many higher conserved currents J^{x) 



(x) - g^C MT and the 



1) 
2) 

3) 
4) 

■5) 



P 



N** q-(6,z) 

Nf^ q+ {e,z) 



N n E e±8t E ^ Q-&Z) 



(2) 



P T+ '(0,z) 



P J -(0,z) 



n 

±Nfr E e ±e ' E e 

i=l i=l 



Lz, 



(L = ±l,±3,...). 



where g± (0, *) " X e * : Y.'" Y. • 

i—1 i—1 i—1 i—1 

The identification with the operators is made by comparing the exact results 
with Feynman graph expansions. Properties as charge, behavior under Lorentz 
transformations etc. will also become obvious. 

The article is organized as follows: In section || we recall some formulae of 
[^| which we need in the following. In section |§| we present a general formula 



3 The symbol J\f refers to norma 



products of local quantum fields. In perturbation theory 



they arc defined by Zimmerman's [ 231 subtraction method, for example 
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for solitonic form factors. In section [| we discuss several explicit examples and 
perform perturbative checks for two- and four-particle form factors. As an ex- 
ample we also investigate the asymptotic behavior when one soliton momentum 
goes to infinity and compare the result with typical bosonic behavior. Section 
^| contains the derivation of a general crossing formula. Again LSZ-assumptions 
are used. Also the charges of the infinitely many higher local conservation laws 
on all states are calculated. In sections [| we formulate the 'bootstrap' prin- 
ciple and introduce the 'bound state intertwiners' (see also Q). Using these 
techniques we derive in section |?] from the pure soliton anti-soliton form factors 
the mixed breather soliton and the pure breather form factors. Section || con- 
tains conclusions and an outlook. Several proofs and explicit calculations are 
delegated to appendices. 



2 Recall of formulae 

In this section we recall some formulae which we shall need in the following 
sections to present our results. All this material can be found in || including 
the original references. Coleman |2f| had shown that the sine-Gordon and the 
massive Thirring model are equivalent on the quantum level. The corresponding 
classical models are defined by their Lagrangian's 

C SG = i^9V+"J(cos^-l) 

C MT = ^(ijd - M> - \g fj» , {f = ?7 M V0 ■ (3) 



2.1 The S-matrix 

The sine-Gordon model alias massive Thirring model describes the interaction of 
several types of particles: solitons, anti-solitons alias fermions and anti-fermions 
and a finite number of charge-less breathers, which may be considered as bound 
states of solitons and anti-solitons. In this work we will concentrate on states 
consisting of solitons and anti-solitons. Integrability of the model implies that 
the n-particle S-matrix factorizes into two particle S-matrices. In particular 
scattering conserves the number of particles and even their momenta. The two 
particle S-matrix contains the following scattering amplitudes: the two-soliton 
amplitude a(9), the forward and backward soliton anti-soliton amplitudes b{6) 
and c(0):| 

, , srahO/v .„, ,„. sinhi7r / V , , 

K0)= . hr ' , o fl , c(9) = —-. ^-o fl , 

smn(i7r — v)/v smh^iTV—vj/v 

f°° dt sinh|(l-^)i . 

a(6) = exp / ^- V" sinhf— . (4) 

Jo t sinh -^vt cosh in 



4 This S-matrix has been obtained first by Zamolodchikov [gfj extrapolating a semiclassical 
result and by means of the 'bootstrap program' using the Yang-Baxter relations in Jl(J. 
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The parameter 8 is the absolute value of the rapidity difference = |#i — 6 2 \ 
where 8i are the rapidities of the particles given by the momenta pi — M sinh 8i . 
The parameter v is related to the sine-Gordon and the massive Thirring model 
coupling constant by 

i? IT 



8tt - I3 2 tt + 2g 

where the second equality is due to Coleman |25| . 

We list some general properties of the two-particle S-matrix. As usual in 
this context we use in the notation 



v 



l...n 



G V 1 - 71 = V! <g> • • • ® K 



for a vector in a tensor product space. The vector components are denoted by 
D ai " ,a ". A linear operator connecting two such spaces with matrix elements 

Aal.'.'.a™' is denoted by 

AxZn ■ V 1 " 71 -> V 1 '-"' 

where we omit the upper indices if they are obvious. All vector spaces Vi are 
isomorphic to a space V whose basis vectors label all kinds of particles (here 
solitons and anti-solitons, i.e. V = C 2 ). An S-matrix such as Sij acts nontrivial 
only on the factors Vi <S> Vj . 

The physical S-matrix in the formulas above is given for positive values of the 
rapidity parameter 8. For later convenience we will also consider an auxiliary 
matrix S (81,62) regarded as a function depending on the individual rapidities 
of both particles 8\ , 82 or some times also on the difference 8\ — 82 

q m ay-q m <n _ J (^Ml^ - 2 |) for 9 1 >9 2 

with a taking into account the statistics of the particles. It is a diagonal matrix 
(T12 with entries —1 if both particles are fermions and +1 otherwise (see J27J). 
The matrix 8(81,82) is an analytic function in terms of both variables 8± and 
82. This follows from unitarity S'S = 1 and the fact that the physical S-matrix 
is the boundary value of a real analytic function S(s + ie) as a function of the 
Mandelstam variable s = (pi + P2) 2 such that S'(s + ie) = S(s — ie) or 

St(0) = S~ 1 (6) = S(-8) (5) 

The auxiliary matrix S12 acts on the factors V\ ® V2 and in addition exchanges 
these factors, e.g. 

512(0) : Vi<g>Vn-+V2®Vi. 



It may be depicted as 



Here and in the following we associate a rapidity variable 0, € C to each space 
Vi which is graphically represented by a line labeled by 9i or simply by i. In 
terms of the auxiliary S-matrix the Yang-Baxter equation has the general form 



S'l2(^12) 513(013) 523(023) — 5*23(023) 5i3(0i3) ^12(^12) 

which graphically simply reads 



(6) 





Unitarity and crossing may be written and depicted as 



£21(021)512(012) = 1 : 



(7) 



12 12 



5i 2 (0i - 2 ) = C 22 5 21 (0 2 + »tt - 0i) C 22 = C 11 S 2l (9 2 - (9 1 - in)) C 11 (8) 






where C 11 and Cjj are charge conjugation matrices. For the sine-Gordon model 
the matrix elements are C Q/3 = C a p = 5 a p where (3 denotes the anti-particle 
of (3. We have introduced the graphical rule that a line changing the "time 
direction" also interchanges particles and anti-particles and changes the rapidity 
as 9 — > 9 ± in. We depict this as 



9 — in 











■ in 



Similar crossing relations will be used below to formulate the properties of form 
factors. 

Finally we note a property of the two-particle S-matrix 



5^(0) = S s a S] 



7 

a u /3 



(9) 



which turns out to be true for all examples. This means that 5 for zero mo- 
mentum difference is equal to minus the permutation operator. 
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2.2 Form factors 



For a state of n particles of kind a.i with momenta pi and a local operator 0(x) 
we define the form factor functions O ai an (6±, . . . , 9 n ) by 

( | 0(x) | ai(pi), ■ • • , «n(Pn) ) m = e-^+-+P^O«(6) , for 0! > • • • > B . 

(10) 

For all other arrangements of the rapidities the functions Og_{§) are given by 
analytic continuation. Note that in general this analytic continuation does not 
provide the physical values of the form factor. These are given for ordered 
rapidities as indicated above and for other orders of course by the statistics of 
the particles. The Og_(9) are considered as the components of a co- vector valued 

function 0i... n (0) € VL.n = (V^-^Vhich may be depicted as 

Ol...n(g) = ( ° ) • 

Now we formulate the main properties of form factors in terms of the func- 
tions 0i...n(#) which follow from general LSZ-assumptions and "maximal ana- 
lyticity". The later condition means that Oi... n (8) is a meromorphic function 
with respect to all 0's and all poles in the 'physical' strips < Im Oij < ir (9ij = 
6i — 9j i < j) are of physical origin, as for example bound state poles as discussed 
in section ^. 

Properties: The co-vector valued function 0i. ..„(#) is meromorphic with re- 
spect to all variables 9\, . . . , 6 n and 

{%) it satisfies the symmetry property under the permutation of both the 
variables 0i, 9j and the spaces i,j at the same time 

0...ij...(- ■ • J 0»> ®ii • ■ •) = O—ji -X- ■ ■ ) ®h ®ii ■ • • ) Sij(9i — 9j) 

for all possible arrangements of the 9's, 

(ii) it satisfies the periodicity property under cyclic permutation of the rapid- 
ity variables and spaces 

Oi,.. n (9i,92, ■ ■ ■ , 6 n , ) = 02...7il(#2, ■ • • , 9 n , 6\ — 2iri)o ox 

(iii) and it has poles determined by one-particle states in each sub-channel. In 
particular the function GgJJP) has a pole at 6*12 = iir such that 

Res 0i...„(0i, . . . , n ) = 2i Cia 3 ... n (9 3 , . . . , 9 n ) (1 - S 2n ■ • ■ S23) 

9 12 =17T 

where C12 is the charge conjugation matrix. 
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(iv) If the model also possesses bound states, the function Oa(9) has additional 
poles. If for instance the particles 1 and 2 form a bound state (12), there 

(12) (12) 

is a pole at On = iu\ 2 ' (0 < u\ 2 ' < n) such that 

Res i2 0ia...„(0i,02,...,0») = O(i2)...„(0(i2),--. A)^lf 2 2) 

012=Mii2 

(12) 

where the bound state intertwiner T\ 2 and the relations of the rapidities 
Qi,02, 9(12) an d the fusion angle u{ 2 ^ will be discussed in section [| below. 

(v) Since we are dealing with relativistic quantum field theories Lorentz co- 
variance in the form 

Oi.„ n (0i +n,...,e n + n) = Oi...„(0i, . . . , e n ) 

holds if the local operator transforms as O — > e s/1 where s is the "spin" 
of 0. 

In the formulae (i) the statistics of the particles is taken into account by 
S which means that S12 = — S12 if both particles are fermions and S12 = S12 
otherwise. In (ii) the statistics of the operator O is taken into account by 
a 01 = — 1 if both the operator O and particle 1 are fermionic and a 01 = 1 
otherwise. 

The properties (?) — (iv) may be depicted as 




As was shown in || the properties (i) — (Hi) follow from general LSZ-assump- 
tions and "maximal analyticity" . The bound state form factor given by (iv) was 
discussed in Q for special cases. In section ^ we investigate the general case 
and show that the bound state form factor is consistent with the 'bootstrap 
principle' which means that it also satisfies (i) — (Hi) if the constituents do. In 
section we derive from the same assumptions a general crossing relation which 
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implies (ii) and (hi). Conversely, it has been shown ^8|, Q that functions 
satisfying the properties (i) — (v) and the general crossing relation represent 
local operators i.e. they are form factors of x-dependent operators 0(x) which 
commute (anti-commute) for space like differences of the arguments. 

We will now provide a constructive and systematic way of how to solve the 
properties (i) — (v) for a co-vector valued function /i... n (£) once the scattering 
matrix is given. These solutions are candidates of form factors. To capture 
the vectorial structure of the form factors we will employ the techniques of the 
algebraic Bethe Ansatz which we briefly explain now. 

2.3 The 'off-shell' Bethe Ansatz co-vectors 

As usual in the context of algebraic Bethe Ansatz we define the monodromy 
matrix as 

2i...n,o(0j #o) = Sio(8i — 9q) 820(82 — 6q) ■■■ S n o(dn — do) (11) 











1 


2 


n 






It is a matrix acting in the tensor product of the "quantum space" V " = 
Vi <£>■■•<£> and the "auxiliary space" Vb (all Vi = C 2 = soliton-anti-soliton 
space). The sub-matrices A,B,C,D with respect to the auxiliary space are 
defined by 



T\... n ,o(S., z ) 



Al.^z) Bi...„(0,2) 
Ci...„(g,js) Dl.^&z) 



A Bethe Ansatz co- vector in V\... n is given by 

*i. ..«(£,!) = ni...nC 1 ... n (6,z 1 )---C 1 ... n (6_,z m ) 



(12) 



where z — (z±, . . . , z m ). Usually one has the restriction 2m < n and the charge 
of the state is q = n — 2m = number of solitons minus number of anti-solitons. 
The solitons are depicted by | or <— and anti-solitons by J. or — >. The co- vector 
is the "pseudo-vacuum" consisting only of solitons (highest weight states) 

fti...n =T ® •• - ® T • 
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It satisfies 



Oi... n 5i... n (0,*) = 

n 

£ll...nM...n(JLz) = II a(6i - z)fLl...n 

i=l 
n 

n 1 ... n D 1 ... n {e,z) = n K^i - *)0i...„ . 

i=l 

The eigenvalues of the matrices A and £), i.e. the functions a = —a and b = —b 
are given by the amplitudes of the scattering matrix In the following we 
use the co-vector ^i...n(&, z) in its 'off-shell' version which means that we do 
not fix the parameters z by means of Bethe Ansatz equations but we integrate 
over the z's M, M, Ejl. 



3 The general form factor formula 

In this section we present our main result. We derive a general formula in terms 
of an integral representation which allows to construct form factors i.e. matrix 
elements of local fields given by eq. (0). More precisely, we construct co- vector 
valued functions which satisfy the properties (i) — (v) on page ||. 

As a candidate of a generalized form factor of a local operator 0(x) we make 
the following Ansatz for the co- vector valued function 



O x ... n {6) = J Ce dzi • ■ • J Ce dz m h{e,z)p°(9,z) *i... n (g,2) 



(13) 

with the Bethe Ansatz state ^i... n (S_,z) defined by eq. (|12|) . For all integration 
variables Zj(j = 1, . . . , m) the integration contours Cg_ consists of several pieces 
(see figure [y): 

a) A line from — co to oo avoiding all poles such that Im 6 L — it — e < Im Zj < 
Im#i — 7T. 

b) Clock wise oriented circles around the poles (of the <j>{0i — Zj)) at Zj = 0, 
(i = l,...,n). 

Let the scalar function (c.f. g) 

n rn 

h(8,z)= n f(%)nn^- 2 i) n ^-^o. ( u ) 

l<z<j<n i=X j=X l<.i<j<m 



be given by 



T(Z) = ^R' (* + »*) (15) 



and 



1 „ f 00 (it sinhi(l-z/)t 1 - coshifl - 0/(i7r)) . , 

F{9) = sin— exp / ^ 4= (16) 

2« Jo * sinh.^i't cosh^t 2smhi 



11 



O u n 



+ in(2v - 1) 2 + in(2v - 1) 



o l + m{2v- 1) 



• t>l + ITT 



n + ITT • V 2 + ITT 

01 + iir(y — 1) 
'„ + iix{v - 1) 6» 2 + i7r(i/ - 1) Q ^ 



O 



O «1 - 17T 



O - »7T o 2 - IK [ 
» ' 

• — 27TI • 02 — 2m 



* 9i — inv 

• 0i - 2ni 



Figure 1: The integration contour Cg (for the repulsive case v>\). The bullets 
belong to poles of the integrand resulting from u{0i — Uj) <p(0i — Uj) and the small 
open circles belong to poles originating from t(0i ~ Uj) and r(0i — Uj). 

The function F(0) is the soliton-soliton form factor. It is a solution of Watson's 
equations 

F{6)=F{-6)a{6)=F{2m-9) (17) 

with a(0) = —a{0) where a(0) is the soliton-soliton scattering amplitude. It is 
the uniquely defined 'minimal' solution || which has no poles and no zeroes in 
the 'physical strip' < Im 6 < ix and at most a simple zero at = 0. 

Remarks: 

• Using Watson's equations ( p7| ) for F(z), crossing (^) and unitarity (^) 
for the sine-Gordon amplitudes one derives the following identities for the 
scalar functions c/)(z) and r(z) from the definitions ( |l5| ) 

1 s o,(z — 2lTl) , . . . 

<f>(z) = 0(*7T -z) = T —(j>{z- *7T) = ^ - r ' 4>{Z ~ 2lTi) , (18) 

b(z) b(z) 

/x , * b(z) a(2iri — z) , „ . 

r (z = r(-z = T7— ( r(z - 2m) 

a{z) byZiri — z) 

where b(z) is the soliton-anti-soliton scattering amplitude related to a(z) 
by crossing b(z) = a(iir — z). 

• The functions <fi{z) and t(z) are of the form 

1 f 00 cftsinh±(l - v)t (coshi(± - z/(wr)) - l) 

4>{z) = const. —— — exp ' 



sinh z Jo t sinh | i/£ sinh t 
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t(z) = const, sinh Z smb. Z / V 



• The function h(6_,z_) and the state ^i... n (0,z) are completely determined 
by the S-matrix. 

In contrast to the functions F{z), 4>(z) and t(z) the 'p- function 

' P°(i,z) in 

the integral representation ( |l3| ) depends on the local operator O(x), in particular 
on the spin, the charge and the statistics. The number of the particles n and the 
number of integrations m are related by q = n — 2m where q is the charge of the 
operator 0(x). The p- functions is an entire function in the Zj (j = 1, . . . , m) 
and in order that the form factor satisfies the properties (i) — (v) it has to satisfy 
the following 



Conditions: The p-function p„(9,z) (where n is the number of particles and 
the number of variables 9) satisfies 

(i 1 ) Pn(0,,z) is symmetric with respect to the #'s and the z's. 

(ii 1 ) p^(0_,z) = 0oiPn(- ■ ■ 2tt«, ■■■,%) and it is a polynomial in e ±Zj (j = 
1, . . . , m). 

The statistics factor oo% is —1 if the operator O(x) and the particle i are 
both fermionic and +1 otherwise. 



(Hi') 



p£(0l = 6 n + in, ~6,6 n ;z.,Z m =6 n ) = - p%_ 2 (9, 1) + pW (9) 

m 

p°(d x = e n + nr,e,e n ;z, z m = e x ) = a ol -p°_ 2 (e,z)+p {2) (0) 

where 0_ = (6>2, ■ ■ ■ , n -i), 1 = (zi, ■ ■ ■ Zm-i) and where p^- 1,2 \9) are inde- 
pendent of the z's and non- vanishing only for charge-less operators O(x). 
The constant x depends on the coupling and is given by (see formula (B.7) 



x=-(F'(0)) 2 /n. (19) 

(iv') the bound state p-functions are investigated in section |^ 
(v 1 ) Pn(0 + H,z + u) = e s, *Pn(0,z_) where s is the 'spin' of the operator 0(x). 

As an extension of theorem 4.1 in Q we prove the following theorem which 
allows to construct generalized form factors. 



Theorem 1 The co-vector valued function Oi...„(9) defined by f u fulfills the 
properties (i), (ii) and (Hi) on page |^ if the functions F(0),<f>(z) and t(z) are 
given by definition filfy - and if the p-function p® (9, z) satisfies the con- 
ditions (i') — (Hi 1 ). 
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Proof. The properties (i) and (ii) follow as in [m. The proof of (iii) is also 
the same as in Q, if the functions p^(0) and p^(9) in (iii') vanish. For the 
case of charge-less operators they are in general non- vanishing. Then the proof 
of (iii) in the form of 

ReS Ol...n(#l, . . . , n ) = — 2i Ci„ 02...n-l(#2, • • ■ , 

X (I2...71-I — S2n • ' 1 5n-lra) 

has to be modified as follows: Considering the z m -integration as in || the terms 
involving Pn-2(]l,z) yield the desired result and those proportional to p^ 1:2 ^(0_) 
yield terms with a factor l2...n-i + $271 1 ■ ' S n —i n which, however, vanish due to 
the following lemma. 

Lemma 2 The integral given by ( j7j| j and vanishes if the p-function p(9_, z) 
is independent of the integration variables Zi and if the number of particles n 
and the number of C-operators m are related by n — 2m which means that the 
charge q = n — 2m vanishes. 

This lemma is proven in appendix 
Remarks: 

• The number of C-operators m depends on the charge q = n — 2m of the 
operator O, e.g. m = (n— 1)/2 for the soliton field tp(x) with charge q = 1 
and m = n/2 for charge-less operators like iptp or the energy momentum 
tensor . 

• Note that other sine-Gordon form factors can be calculated from the gen- 
eral formula ( [l3|) using the bound state formula (iv). 

• The general representation of form factors by formula ( |l3|) is not specific 
to the sine-Gordon model. It may be applied to all integrable quantum 
field theoretic model. So the main task is to solve the corresponding Bethe 
Ansatz. 



4 Examples 

In this section we propose the p-functions corresponding to some charge-less lo- 
cal operators. For two-particle form factors the single integration is performed 
explicitly and compared with known results of || and with Feynman graph 
expansion in lowest order. For 4-particle form factors the double integrals are 
calculated approximately by expansion in small couplings to lowest order. The 
results are again checked against Feynman graph expansions. In all cases agree- 
ment is obtained. 
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For the massive Thirring model with the Lagrangians and field equation 
C MT = ^(ijd - M)i> - \gfi» , (f = ~tf^) 

we are looking for all matrix elements of the quantum operators corresponding 
to the following classical fields 

1) tp(x)tp(x) 

2) ^(x)j 5 i(j(x) 

3) j^(x) = x/j(x)j^ip(x) the topological (electro magnetic) current 

4) T IJ,u (x) = |t/i7 m d v tp — g^ v L MT the energy momentum tensor. 
The light cone components of this tensor are 

T ±± = T oo ± 2T oi + T n = ^ 7 ±lgi^ 

T+- = T-+ = T m - T 11 = i> 1 +l -8^^ - gfjn = MW 

where d = d° ± d 1 and 7 ± = 7 ± 7 1 . For the last equality the field 
equation has been used . 



5) The higher conserved currents |32j, |33j for L = 3, 5, . 



(d+) L 4>i + h.c. + o(^) / i>i 



A second set of higher conserved currents is obtained by exchanging d + <-> 
d~ and ip\ «-> tp2 which we associate to L = —3, —5, .... 

4.1 Examples of 'p-functions' 

In this subsection we propose the p-functions for various local operators. Since 
the charge of the operators which we consider is zero, the number of integrations 
m and the number of particles n are related by m = n/2 and form factors are 
non-vanishing only for even number of particles n — 2,4,.... We consider 
p-functions of the form 



p?(g,2) - n° [ P °(pnJ2 eLZj +p°(pnJ2^ LZ3 I ( 2 °) 

where P M is the total energy momentum vector of all particles. The integrals 
in ( pr] ) converge for L < + l)(n/2 — m + 1) + 1/v. For large values of L 
the form factors are defined in general as analytic continuations of the integral 
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representation from sufficiently small values of v to other values. Obviously the 
p-functions ( p0[ ) satisfy the conditions (i 1 ) — (Hi') on page |l3|. From the property 
(in') we obtain the recursion relation for the normalization constants 

«?=ec => ^=JV?i* m " 1 . (21) 



The absolute normalizations follow from the two-particle form factors (see eq. (|2 
below). 

We propose that the p-functions of equations (0) on page |^ are associated 
to the local operators Af [Ipip] (x), Af [r^y 5 ^] (x),p 1 (x), T ±± (x), T+~(x) and 
Jl(x) where ± denote the light cone components (e.g. — j°±j l ). The vector 
operator j^(x) = Af [ ■07 AI '0] (x) is the topological (electro-magnetic) current, 
T^ v (x) is the energy momentum tensor and the J^(x) are the higher conserved 
currents. The normalization constants are obtained in the next subsection, 
where we calculate the exact two-particle form factors. The fundamental sine- 
Gordon bose field <p(x) which correspond to the lowest breather is related to 
the current by Coleman's formula 

e^d uV = -jf or d ± <p = ±jj ± . (22) 

This implies the following representation for the p-function 

_ . / n n \ — ^ / n m n m \ 

x&z) = n^ £e«£ e -« E^E eZ +E^E^ ■ 

^ \i=i i=i / \i=i i=i i=i i=i / 

Using the integral representation ( [j"3| ) with these p-functions we calculate 
in the following subsections the exact two-particle form factors and the four- 
particle form factors in lowest order with respect to the coupling g. 

Conservation of higher charges The higher currents satisfy in terms of 
matrix elements Ox... n (§) = ( | O \ pi, . . . ,p n )™ n for all L G Z the equation 

d + J^(x)+d-J+(x) = 0, 

such that the higher charges are conserved 



Proof. ^,From the definition we get the correspondence of operators and p- 
functions 

(n n n n \ m 

E^E e ~ e! -E e ~ e! E ee! E^ = °- 
»=i i=i i=i i=i / i=i 

■ 

In the following subsection we calculate the charges on 1-particles states 
and in section ^| on arbitrary n-particles states. It turns out that for even L 
the charges vanish as in the classical case. The energy momentum tensor T^ v 
is given by and the momentum operator by P° ± P 1 = Q±\. 
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4.2 Examples of two particle form factors 

Two particle form factors may in general be obtained by diagonalization of the 
two-particle S-matrix || . For several examples we shall show of p- functions that 
the integral representation in this case reduces to known results for two particle 
form factors for specific operators. This allows us to confirm the association 
of p-functions and local operators as proposed above. Also the normalization 
constants can be calculated. Moreover we check our results in lowest order of 
perturbation theory. 

We consider the form factor for charge-less operators 

O 12 (6 l7 9 2 ) = (0\O\ Pl , P2 )? 2 . 

Non-vanishing matrix elements contain one soliton and one anti-soliton. Choos- 
ing n = 2 and to = 1 in the general formula ( |l3| ) we obtain 

Ou{6) = F(9 12 ) f dz]l<f>(e i -z)p (6,z)n 12 C 12 (e,z). (23) 
Jet i=l 

The integration can be performed. As a special case we take a p-function of the 
form (J2G) for L = 1 and prove the following lemma: 



Lemma 3 For the simple p-functions e ±z the integral in ( pq j can be performed 
yielding the result 

f± (e) = 2siDh ^i2 e ±i(e 1+e2 ) ( ± U(0i2)E+ _ f-(e 12 )E^ \ 
12{ -' vk \ coshi6» 12 sinh±6> 12 / 

where £ ,± = (s ® s ± s ® s) is the symmetric (anti-symmetric) soliton anti- 
soliton state and the constant x is defined in eq. (fH|). The functions f±(&) 
are the positive and negative C-parity two-particle sine- Gordon form factors 
calculated in j^/ 

(/ (0) / (0)) f ta " nh l( i7r ~ 9 ) 1 \ F(9) 

' \ sinh ^- (^TT — #) ' cosh ^- (in — 0) ) 



with /i? (8) given by eq. (jii 
Proof. We consider the expression 

I* = ^4/ dzI(z)e±*nC(9,z) 

C[0\2) J Co. 

°JM I ' t w <*. { m - .) ^ - .) + 4^4 , 

[a(Ui — z) a{0 2 — z) a(a\ — z) 



Co 



with i{ z ) = rii=i,2^(^ - z ) . = wmq) and 

C(§-: z) — C(^7 z ) n<=i 2 l/ a (^i — z )- Inserting the identities (which follow from 
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Yang-Baxter relations for the soliton S-matrix) 

a(9 12 ) a(6i - z) a(9 2 - z) b(9 1 - z - 2m) b(9 2 - z) 



c(9 12 ) 



c(6i - z) c(9 2 - z) c(6i - z - 2m) c(9 2 - z) 
a(0i - z) a(9 2 - z + 2m) b(0 1 - z) b(8 2 - z) 



c(6»i - z) c(9 2 -z+ 2m) c(9 x - z) c(8 2 - z) 
into the two components of the integral, respectively and using the shift property 



z — 2ni) 



b(0-z) - 



we obtain 



/-/ 


dzl(z)e ±z | 







Ic(0!-Z) 



*) 



a{9 2 - z + 2m) 
c(6 2 - z + 2m) 



= -2m 



I(z)e 



±z 



h - z) a(8 2 - z + 2m) 

■S®S h — -f 80S 



c(0i - z) 



c(9 2 - z + 2m) 



There are no poles inside the integration contour. However, there are contribu- 
tions at ±oo. With the asymptotic formulae for Rez — > ±oo 



a{9 - z)/b{9- z) 
b{9 - z) / c{9 - z) 

4,(9 - z) 



e ±tir(l/«-l) 

±1 
2i sin(7r/i/) 
4 



e ±fa(l/i>-l) e ±i(l/«+l)(e-a-t7r/2) > 



which are derived in appendix |b] and with c(9) / a(9) = i sin — / sinh ^ (in — 9) we 
obtain the claim. H 



4.2.1 The exact form factors 

We use the following conventions for the 7-matrices and the spinors 



7° = 



1 

1 



-1 1.7" = -"-,' 



e -e/2 \ / e -e/2 

u( P ) = VM ( e/2 , v( P ) =VMi[ g/2 



For the examples 1) - 4) above we calculate the two-particle form factors 
O sS (6 1 ,6 2 ) = (0\O(0)\p 1 ,p 2 )fs 
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applying lemma || and using the p-functions (^) explicitly]^] 

1) (0\Ar\W>](0)\ Pl ,p2)% = v{0i)u{6 x )f + {e 12 )/v 

2) (OljV^T 5 ^] (0)|pi,Pa)S3 = 0(W«(W-(0i2)/f 

3) (01^(0) |pi,^)*5 = ^2)7^1) /-(M 

4) <o|^(o)b 1)P2 )™ = e(0 a )yu(0i)§(tf-j£)/+(0i 2 )/i/. 

(24) 

For the higher conserved currents see the next paragraph. The normalization 
constants in the expressions for the two-particle form factors of (0) have 
been determined by the following normalization conditions 

1) s (p\M [V^] (0)\p) s = u(6)u(9) = 2M 

3) s (p\Af [^y^] (0)\p) s = u(6)^u(6) = 2p» 

4) s (p\T^{Q)\ P yj' = u{9)^u{9)p v = 2p^p v 

which are the free field values and which are natural due to the corresponding 
'charges' of the operators. The crossing relations and /+(z7r) = v and f-(in) = 1 
have been used. Since there is no 'charge' for AT [ t/j^tp] ( x ) we take N^ 1 * = 

—N$ which implies the natural relation s ( p' \Af ^tjjj 5 t/j] (0) \ p) s ~ u(9')j 5 u(9) 
for 9' k for small couplings. This normalization is also consistent with the 
desired identification (see |25| ) 

a^(i±7 5 h (-)=-^( i -Q- ± ^ ): 

where : • • • : means normal ordering with respect to the physical vacuum (see 
also J34|, [55]). The constant on the right hand side is obtained from the trace 
of the energy momentum tensor calculated in [|4| E35j. For the cases 3) of the 
topological (electro-magnetic) current and 4) the energy momentum tensor the 
given normalization are equivalent to the eigenvalue relations 

dxf{x)\p) s = \p) s 

dxT»°(x)\p) s = p»\ P ) s . 

Finally we obtain from the recursion relation (|2l]) and lemma || the normaliza- 
tion constants: 

N** = -ATfT 5 ^ = ^-W = -^pN? = -iM^—x m (25) 
2v " M 4m! 



5 For the cases 3) and 4) these results agree with those of |2| and [W, respectively, which 
have been obtained by solving the scalar Watson's equations aue to diagonalization of the 
S-matrix. 
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where the constant >r is defined in eq. (pit). Note that these normalizations 
together with the proposed p-functions (0) imply in particular the quantum 
version of the classical field operator relation 

trT(x) = T+-(x) = MM \jjrip] (x) 

for all matrix elements. For the sine-Gordon field Coleman's relation yields 



Vs ,(0u9 2 ) = (O\ V (O)\ Pl ,p 2 )T s = 2 " 1 ! "'- :) 



/3 cosh i 812 



The higher conserved currents: We will use a generalized version of lemma 
^ where e ±z is replaced by e Lz . The two particle form factors then turn out to 
be 

[Jt]n (<Ma) = <°l J t |Pi,Pa>$ = ±N^ sinh|0 12 (E e± ° 4 ) e l L &+°>) 

x k L (0 12 ) ( S f ^ f + (8 12 )E+ , \ /-(M^iO • (26) 
\coshiy12 smhie'12 / 

The asymptotic behavior of <p(0) yields the functions (see appendix |b|) 
k L (8) = cmst^^Si+w^AiAjeft'-M 

i=0 j=0 

where the constants Ai are also given in appendix |b| Only the first term pro- 
portional to f + (0\ 2 ) contributes to the higher charges Ql = J dxJ®. With a 
suitable normalization (N^ k^^Tr) v = M L+1 ) and for odd L = ±1, ±3, 
one obtains 

s(p2\Ql\pi)s = 2nS(pl -p\)\{ [Jt] sg (0i,0i - «0 + [JI] sg (9i,9i - Mr)) 
= ( P2 \ Pl ) (p+) L . (27) 

Note that k^i-K) = for L even. 



4.2.2 Feynman graph expansion 

From the Lagrangian (||) we get the Feynman rules of figure ^. The two-particle 
form factors for charge-less operators of the form 0(x) — AfipTip are given in 
lowest order perturbation theory by the Feynman graph depicted in figure ^. 
We consider several examples for V: 
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k' 1 



■yk - M 



Figure 2: The Feynman rules for the massive Thirring model. 




Pi 



I>2 



Figure 3: The Feynman graph for two-particle form factors of charge- less oper- 
ators. 



1) (0|^|pi,P2>s. = -w(0i)u(0 2 ) = -2Misinh±0i 2 

2) (0\^j 1 5 i;\p 1 ,p 2 } ss = -v{6 1 ) 1 b u(6 2 ) = 2MicosYv\6 12 

3) (0\H k ^>\Pi,Pa)s. = ~ v(e 1 ) 1 ± u(e 2 ) = T 2Mie ± ^+ e ^ 

4) (o|^7+^a+v>|pi,P2>s. = ~v(e 1 ) 1 +u{e 2 )\{ Pl -p. 2 )+ 

= 2M 2 ismh±0 12 e ±( - 9l+e ^ 

Note that the energy momentum tensor light cone components are T ±=t = 
|jV ifr/^d^ip an( l T ±T = MA/" fjpi/A ■ For all these examples we have agree- 
ment with the exact expressions ( p4| ) when g — > or v — > 1. 

4.3 Examples of 4-particle form factors 
4.3.1 Expansion of the exact formula 

We investigate the integral (|l3|) (for n = 4 and m = 2) 

with the scalar function 

4 2 

l<i<3<4 i=lj=l 



<fei / <feaM0,*)p u (£,z)*s* M (0,*) 

Co JC e 
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and the Bethe Ansatz state component 



= (flC(9, Zl )C(6, z 2 )) ssss = H H a & - z j) 



i=l 3 = 1 



x (5(0! - Zl )c{0 2 - z 2 ) + 6(0! - zi)c{0 2 - Z!)c(0i - z 2 )~b{6 2 - z 2 j) . 

with b = b/a, c = c/a for small couplings. We consider first the simple p- 
functions 53j=i e±Zi an< ^ using the Z\ «-» Z2 symmetry we calculate the following 
integral in lowest order with respect to the coupling constant g 



J± = 



/ rfz! / dz 2 (nn tv* - ^ ) ^ - zi ) ^ - **) 

JCf, JCa \i=l 3=1 ' 



x (l + b(9 1 - z 2 )b(9 2 - Zl j) t{z! ~ z 2 ) (e ±Zl + e ±22 ) 



8ig sinh lfl 12 sinh lfl 34 e **(*+fr) * ±6( 2 
n^^inhi^cosh^,)^^^)^ 6 +U W>- W 



The derivation of this result is quite involved. A sketch of the calculation is 
delegated to appendix |c[ 

Finally we use the p-functions (^) and the normalization constants given by 
eq. ( p5| ) and F(9) = —i sinh ^9 + 0{g) and obtain the four particle form factors 



\ in 

I ssss 



O SSss {6) = (0|O|pi,p 2 ,P3,P4; 
for the various operators in lowest order in the coupling g as 

i — i 1 sinh h #10 sinh h #34 , , , v-^ a v- 

\W>] - ssss (0) = -?9M 2 - 2 cosh § <?13 + M E e*> 2 e- 1 

r — k n , x 1 sinh i sinh 4 #34 , , , ^ — ^ „ ^ — ^ fl 

GfrV] _ (g) = ~ 2 9M n 2 34 sinh I (# 13 + M E e* E «"* 

J. Xl< 7 Z J 7 — 1 1 



U<j owoli 2"W i=l i=l 

sinh i#i2 sinh i# 34 . 
U<j cosh 5°v 



11 / 4 

1 sinh -#i2 sinh ^#34 , „ ^ / ±e, 



s - s - ss (g) = TgM ° 2 "° lfl 2 34 sinh \ (9 13 + 9 2i ) £ ^ 

2 ' rUj-coshi^ v . =1 

1 , r2 sinh |#i2 sinh i# 34 1 ! 

ll f<j cosh±0, 3 ._ i ._ i 



[^1 ^ (fi) = cosh I (# 13 + M E « 



e 



[t+-] s _ sss m = -t gM * i-i^iii cosh i (# 13 + # 24 ) e^E' 

2 rii<jCOsh2#y ^ ^ 

which agrees with the Feynman graph result calculated in the next subsection. 
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Pi P2 P3 Pi Pi P2 P3 Pi 

(a) (b) 



Figure 4: Feynman graphs for four-particle form factors of charge-less operators. 
4.3.2 Feynman graph expansion 

The graphs of figure |] give the lowest order contributions to the matrix element 

(0\lpTlp\p 1 ,p2,P3,P4Ys r sSs = — *ifl r <0 | <0r-0 J d 2 xfj^\p 1 ,P2,P3 i Pi)ssSs 

= -ig [(G a + G b ) - (1 «-> 3)] - (2 <-> 4) 

up to terms of order 0(g 2 ). We obtain 

i 

1VP2 +P3 +P4) ~ M 
■ 1 - r l(P2+P3+P4) + M - a - 
2M (P2 +P3+ Pi) - M 

1 

G b = Vi^U 2 V3l^—, : ^7 ru 4 

1{-Pl - P2 - P3) - M 



■ 1 - a - - j(-Pi ~ Pi ~ Pa) + M 

= - l T7VJ v iT u 2 v 3 "f^v 2 v 2 -, ■ ■ ^ —^Tu^ + symm 

2M (pi +p2 +P3) — 



with Ui — u{9i) etc. The relation uu — vv — 2M1 has been used. The terms 
symm 



vanish after anti-symmetrization. We obtain for T = 1,7 5 , 7 ± , 7 P |'<9 5> 



G a = AfitaDh|0 24 - 2 12 -(2<->4) 

COSh i& 3 4 

, sinh ^#14 
G b = Mitanh±0 13 -f-ii - (1 <-> 3) 



Gl = -Mi tanh ±(9 24 nr 1 - (2 <-> 4) 



Gl = -Mztanhffl 13 ' " - (l <-> 3), 
cosh i6» 23 



sinh |6 


>12 


cosh if 


?34 


sinh ~6 


>14 


cosh |< 


'23 


cosh i 


#12 


4 cosh | 




cosh i 


6\4 
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±i(e 1+ e 2 ) 

Gt = ±Mi tanh \6 2i — =- (2 <-> 4) 

COSH i 34 

±i(ei+e 4 ) 

Gf = ±Mi tanh ±8 13 - T - (1 «-► 3), 

COSh 



G» a ° = pM\^U 2 ^ 6 ^ (-a Slnh \° 12 e <4^+^ + e'ifr+o*)) 

\ cosh i6» 34 / 

-(2-4) 

G p ° = P M 2 itanh±9 13 e^+ e *'> (-a — f 14 e g *( 9l+fl «> - e^+SaA 

V COsh ^6*23 / 



which gives after anti-symmetrization the same expression as the one calculated 
from the integral representation of the exact form factors (exchanging 2 <-> 3 
and the sign due fermi statistics). In particular for the energy momentum tensor 

we obtain with T ++ — ipj + tp in lowest order in the coupling g 



[0\T ++ \pi,p 2 ,P3,P4 



\ i n 



_ 1 2 sinh ^6*i3 sinh 1 6> 2 4 cosh |(6>i 2 + 6> 34 ) / g \ 

Furthermore we have agreement with the classical relation for the trace of the 

energy momentum tensor: T H = ip^ ' \d~ — g j^j^ = Miptp since the Feyn- 
man graph calculation implies 

(0\fy + ^r4-MW\pi,P2,P3,P4)sis. = 16 5 A/ 2 sinh i# 13 sinh ^ 24 

= ff(0| 3*3li\Pl,P2,PZ,Pi)Tsss 

in the coupling g as it should be. 
4.3.3 Asymptotic behavior 

We are interested in the asymptotic behavior of form factors when one or more 
rapidities tend to infinity. In perturbation theory for pure bosonic models one 
may use Weinberg's power counting theorem for Feynman graphs j35|[| For the 
exponentials of the boson field O = Me llv this yields in particular the asymp- 
totic behavior 

0»(0i,0a, ■ • • ) = Ox (ft) O n -i(6 2 , ■■■) + 0(e~ Rcei ) (29) 

6 This type of arguments has also been used in |], |36[ [sv], H . 
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as Re 61 — ► oo in any order of perturbation theory. This behavior is also assumed 
to hold for the exact form factors and it was used e.g. in p4j to obtain the 
normalization of exponentials of fields. For fermionic models the asymptotic 
behavior is more complicated. As an example we investigate a component of 
the four-particle form factor of the operators O = ip (l ± 7 s ) tp for the massive 
Thirring model 

0% ss (0) = (0\O ± (0) \p )f Sss = f dz x l dz 2 h{9, z) p ± {9, z) * S5m {9, z) 



with the p-functions of eqs. (|2j|2£]) 

P ± (0,z) = T NY / e ±e *Y, eTZ *- 
After some calculation we finally obtain the asymptotic behavior for Re B\ — ► 00 
Of Ss M) » const* e^- 1 '"^ x f dzh(6',z)q ± {9',z)*s SS (9',z) 

with & = {62,63,64) and 

const^ = const / dz e^^ 1/u±1)z (j){~z) c{-z) 



'Co 

>l h {6!,z) = e-KV^^ei^+D^j e_ 

i=2 



* Sss (0',z) = c(0 2 - z) a{6 3 - z) a{9 4 - z) 

Obviously this is not of the form given by relation (^9|) , in particular the func- 
tions q {6', z) are not valid p-functions since they do not satisfy the conditions 
on page O. This means that they do not correspond to local operators. 



5 "Crossing" 

5.1 General crossing relations 

In order to obtain the general matrix element OM *(0' | 0{x) \ 4>) m of a local 
operator from that case where the state <j)' is the vacuum one uses 'crossing'. 
This means one shifts in the matrix element particles from the right hand side 
to the left hand side. It is well known from the theory of Feynman graphs and 
more general also from LSZ-reduction formulas, that these shifts are related 
to analytic continuation. We will now derive by means of LSZ-assumptions 
and 'maximal analyticity' a formula^ which gives a general matrix element of 

7 In y a similar formula was proposed which differs from.the results of this paper by sign 
factors. Our proof of the crossing formula follows that of te4|. 



25 



a local operator 0(x) in terms of an analytic continuation of the form factor 
function 0\„, n (ff). As a generalization of the co- vector valued function 0i. ..„(#) 
we introduce the short notation 0/(0j;0j) given as follows. 

Let the array of indices / = (ii, . . . ,4m) denote the factors of the tensor 
product of vector spaces Vj = Vi 1 ,,,i, I , = Vi ± ® • • • ® and correspondingly for 
J. For a local operator 0(x) and for ordered sets of rapidities il > ■■■ > 9i m 
and 6'^ < ■ ■ ■ < 0j m we write 

Ojie'y,^) := o « t 0| J |(^ |j| ),...,j 1 (^ 1 )|O(0)| ll (p il ),...,z m fe |I| )r (30) 

where 0j = (8 h , . . . , 0; m ) and Sj = 0$, , 7| )■ The function 0/(0^; 0, ) 

intertwines the spaces V/ — ► Vj and may be depicted as in figure 0. Similar to 



h 3\j\ 




Figure 5: The general matrix element of a local operator. 

C>i. ..„(#) this function is given for general order of the rapidities by the symmetry 
property (i) for both the in- and oi/i-states which takes the general form: 

For the connected contributions this follows again from analytic continuation 
#ii > 6i 2 ~^ &ii < @i2 which can be proven similarly as for 0i. ..„(#) ||. For the 
disconnected contributions this may be considered as a convenient definition. 
As a generalization of the two-particle S-matrix (including the statistics factor) 
Si 2 (01,02) we have introduced the more general object Sj(0j) given by the 
following definition. 

Definition 1 Let J = tt(I) be a permutation of I. Then Sj(0_j) is the matrix 
representation of the permutation group <Sm generated by the simple transposi- 
tions o~ij : i <-> j for any pair of nearest neighbor indices i,j £ I a^ 

a ij ^S^< I \e i ) = S ij (9 ij ) 

Because of the Yang-Baxter relation and unitarity of the S-matrix the represen- 
tation Sj^(6j) for all ir 6 5m is well defined. We will also use the notation 

sf M {ej) = s^iOj) 

8 Note that this definition is quite analogous to that of representations of the braid group 
by means of spectral parameter independent R-matrices. 
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if 7r is that permutation which reorders the array I such that it coincides with 
the combined arrays of K and M which means that 

7r(J) = KM = (fci, . . .,k\ K \,mi, . . . ,m\ M \) . 

As an example consider the case / = (1, 2, 3,4), K = (2, 3) and M — (1, 4) 



'S'l234 (^1; 6*2, &3, #4) 

2 3 14 



<5'l3(^13)5'l2(^12) 

2 3 14 



5' 



1 2 3 



1 



3 4 



Similarly S'[ M (Q_ L ,0_ M ) is denned as an inverse by the formula 
S'[ M (d L ,8 M )Sj M (9j) = 8 j where 5j is the unit matrix in the vector space Vj. 

Analogously to eq. ( |30| ) for the general matrix Oj(8'j;8j) for local operators 
0{x) we write for the unit operator 

|A 1 = ™ t (n|iv|(K w ) ) ---^i(^J|m 1 (p ro J,...,m| M |(p m , M| )r 

I M I 



if the rapidities are ordered as 8 mi > ■ ■ ■ > 6 m]M ^ and 8' ni < ■ ■ ■ < 0' and if 
N is the completely reordered array M. This object has obviously no analytic 
properties, since it is completely disconnected. However, we define it for other 
orders of the rapidities again by the form factor property (i) 

This implies in particular that 

\M\ 



> 



> 



i=l 

9 m | M l and 9' ni > ■ ■ ■ > 8' n Here is the unit matrix in the 



for 

space Vm — Vn ■ For example if 8\ > 6 2 and 8[ > 8' 2 one has the two cases 



\ a 'P' la* 



L af3 



) = m { 



m {f3 l {d' 2 ),a l {8' 1 )\a(8 l ),P{d 2 )Y 

2 

Saa' 8/3/3' 



i=l 



{a'(8' 1 ),p'(8' 2 )\a(8 l ),p{8 2 )Y 

2 
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Theorem 4 For any local operator 0(x) the 'intertwiner valued' function 
&j (Jt'ji (ti) defined by eq. pY^ satisfies the general crossing relations^ 



(31) 

LuJV=J 
KUM=I 

LUJV=J 
KUM=I 

where K, L, M, N and K ,6 L ,6 M ,6 N are defined analogously to I andOj. How- 
ever, L = ■ ■ ■ , h) and 8^ — , ■ ■ • ■> ^fj where the bar denotes the anti- 
particles and C LL is a multi-particle charge conjugation matrix. The general 
crossing relations are depicted in figure ra. 




Figure 6: The general crossing relations 

Proof. In || the crossing formula was proven for the case of only one out-going 
particle | J\ — 1 using LSZ-reduction formulae and the assumption of maximal 
analyticity. General LSZ-reduction formulae for bosons take the form 

out ((j)' | O \ p,^Ya = ° ut ((t)'\a a ut \p)0\(t)) m 

+ i J d 2 x out (cb' I T [OjUx)] | (f>) m e~ lpx 

9 These crossing formulae are generalizations of those in |24]. 



28 



ou < i ((j)',p'\0\<p) m = mt (<f>'\Oa™(p')\ct>) m 

+ t J d 2 x out (<f>'\T[Oji(x)} \. 

and for fermions the form 

"* (<f/\0\p,<l>) t = a 0a out W | a™* Hp)O\0) m 

j2 „ out i it 



+ ij d 2 x ° ut ( </>' | T [Oj a (x)] | 0) m u(9) e~*P x 

°t{4>\p' | O | 4>) m = uo & \Oa™(p') | 4>) m 

-ia a& J d 2 x out (<f>' | T [Oj a (x)} \cj)) m v{9') e^' x . 

Using these relations one can derive similar crossing formulae for the general 
case that 4>' ^ 0. In these more general formulae one particle (from several 
particles) in the owi-state is crossed to the in-state. Iterating these relations 
one obtains the general crossing relations. H 

Remarks: 

1. Note that the equivalence of both formulations of the crossing relation 
follows from the properties (i) — {Hi). 

2. As was shown in || the properties (ii) and (Hi) on page [l3| follow from 
the crossing formula for \J\ = 1. 

5.2 Eigenvalues of higher charges 

In this subsection and as a simple application of the general crossing formula 
we investigate the higher conservation laws given by the currents J^(x). 

Proposition 5 In terms of the matrix elements 
Oi... n (0) = ( | O . . . ,p n )i", n the higher charges 

Ql= J dxJl(x) = J dx\ (J£(x) + J^{x)) 

satisfy for odd L = ±1, ±3, . . . the eigenvalue equation 



QL-^2(pt) L )\pi,...,Pn)t.n = 



for suitable normalizations constants N^ L . For even L the charges vanish as in 
the classical case. 
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Proof. As a generalization of equation for n = 1 we prove for arbitrary 
n = \I\ and n' = \J\ that 

n 

lQL]i(e';9)=Y,{pt) L i J i(&';6)- (32) 

i=l 

First we show that for ?i + n' > 2 the connected part of the matrix element 
[Ql}\ "n (§L>§) vanishes. This connected part is obtained by the analytic con- 
tinuation [j£] , , (6L+ iK,9). From the correspondence of operator and 
p- function 



« / n n \ t 

/ dxJl (x) ~ 2n6(P> - P)NX n Ef E ^ ~ E e±Sl E ^ 

J ' i=l i=l / i=l 



the claim follows since for n' + n > 2 there are no poles which may cancel the 
zero at P' = P where pM — Y^Pi ■ Note that only for n = n' = 1 the factor 
1/ coshi6'i2 in ( ]26| ) cancels the zero. Therefore contributions to (32) come from 
disconnected parts which contain (analytically continuated) two-particle form 
factors (c.f. eq. (p7|)) 

IQiM {e j ,e i ) = {pi) L m j A)- 

It follows that in the general crossing formula only those terms with K — 
{i} , M = I\{i},L = {j},N = J\ {j} contribute 

n n 
i=l 3=1 



Y,(pt) L i J Am 



i=l 

where the obvious relation £"=i S[ N (9! L ,§! N ) l^fl,-,^) l^(^^Af) 5f M (£) 
1/ nas been used. ■ 



6 Bound states 

Before we define and investigate the properties of bound state form factors 
we recall some facts on bound state S-matrices and define the "bound state 
intertwiners" . The two-particle S-matrix satisfies real analyticity ([5J) unitarity 
(0), crossing (||), the Yang-Baxter relations (||) and the permutation property 
at vanishing argument (||). 

Let the two particles labeled by a and (3 of mass m a and mp, respectively 
form a bound state labeled by 7 of mass m 7 . The mass of the bound state 7 is 
given by 

TO 7 = Jm 2 a + nip + 2m Q m (3 cosu^ /3 , (0 < u^p < 7r). 
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where v? aa is the so called fusion angle and iu^g is equal to the pure imaginary 
relative rapidity of the constituents a and (3. The two-particle S-matrix may be 
diagonalized 

e 

where the projections onto the eigenspaces (labeled by e) are given by the in- 
tertwiners fp{9) and f e a g{9) with 

e a/3 

The eigenvalue of the S-matrix S(a, [3, 7, 9) which correspond to a bound state 
(a/3) = 7 has a pole at 9 = iu^ g , a fact which will be used to define the 'bound 
state intertwiners'. 



6.1 Bound state intertwiners 



Following the investigations of j27) (see also ]24| ) we use in addition to the 
intertwines y>^g which are defined for all eigenstates e of the S-matrix also similar 
ones r^g(ttXg) which are defined for all fusion angles. They are therefore only 
defined for an eigenstate 7 which correspond to bound states i.e. an to eigenvalue 
of the two-particle S-matrix S(9) which has a pole at 9 = iu^ g . 

Definition 2 The matrix elements ^ 1 a g{9 1 aB ) of the bound state intertwiner 

are defined by the residue of the S-matrix 



i Res S^\e)=T?«\ul )TlJul ) 



i Res 



(3' a' 



a (3 



/3' a' 

U 



where the dual intertwiner is defined by the crossing relation 

13 a 




(34) 



(35) 



with the charge conjugation matrix C (e.g. for the sine-Gordon model C a 
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Remarks: 



1. For the bound state intertwiner given by the matrix elements 



and defined by eq. (E 



we will also use the notation 
(12) 

r& 9) («!?>) = 



where u 



(12) 
12 



itlg. It intertwines the spaces V\ ® Vi and ~V{yi) 



r< 2 2) : vi ® v 2 - v (12) 



(36) 



("12) f 12 ) 

2. Note that the bound state intertwiners T 12 (u 12 ) and the dual ones 



p21 
i (12) 



( u i2^'') are defined for m^ 2 2 ' > 0. In addition one may define the 



'inverse' ones T 2 ^ 2 ^ (u^ 2 " 1 ) with u 2 Y' = —u^' < such that 
r 12 (,, (12 hr (12) („ (12) ) - p fi9i 

i (12)l M 21 J 1 12 ( M 12 ) — -n.2{U) 

7-1(12)/ (12)x-pl2 / (12)n x 
r i2 «2 ) r (12)( U 21 ) = S (12) 

where Pi 2 (12) projects onto that subspace of V\ ® Vq, defined by ( ^6|) on 
is nonzero and 6(12) is the un it matrix in Vm) • The dual 



(12) 



,(12) 



(37) 



which r 



(12)c,.(12) 

12 ("12 



'inverse' intertwiner is again defined by a crossing relation analogously to 



(35) 



3. Obviously the bound state intertwiners are defined only up to some phase 
factors. From eqs. ( ^3|p4| ) follows that the components are given as 



T lp = e ( a >P> 7)* 



Res S{a,B, 7, ( 



1/2 



7 



Here S{a, j3, 7, 6*) is the eigenvalue of the S-matrix Sfg" which corre- 
spond to the bound state 7 and the e's are phase factors. In |^7| it 
was shown that i Res Sj (8) is real where the sign is related to the parity 
of the particles involved. The phase factors e fulfill e(a, f3, j)s(a, /3, 7) 

= — sgn (j, Res S(a, (3, 7, ( 



Example 1 For the sine-Gordon model we fix the phase factors by e(s, s, k) = 
(— l) k and Lp^g = (— l) k <fg s = 1/V2 0. The two-particle S-matrix has poles at 

10 This choice is motivated by th e fa c t, th at with this convention the breather matrix elements 
(0| : <fi k : {0)\b k ) are positive (see Q gj). 
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6 = ivP*> = — kv) which correspond to the soliton anti-soliton bound states 
alias breathers bk- Using the short notation T^ s — T b s g(u^) we obtain 



i* = (-i) fc rt = (-if* 



\ n Res S±(B) 



0=iTr(l-kv) 



1/2 



(38) 



where + and — correspond to even or odd k, respectively. The residues have 
been calculated in \F\j 

Res_ S±(6) = n Res b(6) = (-l) k _ Res c(0)(-l) k 



9=in(l-kv) 



9=in(l-ku) 



—i7r(l — kiy) 
fc-1 

= 2i(-l) fe cotffo/ JJcot 2 ^!/. (39) 



1=1 



All bound state intertwiners T involving solitons are uniquely given by the cross- 
ing relations 



1 s — 1 fcs — 1 fc — 1 sfc — 1 s — 1 ss 



rf = r* fe5 = rf = r s sk = rf = rf = (-i)T* 



Proposition wl on page L?q and the general relation (vA) imply up to a sign 



Res b{6) 

9=m{l-kv) 



Res S ks (6) 

9=t7r^(l+fcr/) 



which may easily checked directly. 



Example 2 For the breather-breather bound states we fix the phase factors 
e(k, I, k+ I) = 1 and (p 1 ^ 1 = 1 (for k + I < 1/u) and get 



kl 



Res S k i{9) 

6=ii:±{k+l)v 



1/2 



where for k < I 



Res S k i(d)=2itimUk + l)v- TT 

9=i7r£(k+i)i/ tan§/ci'- L - L 



tan £Zz/ i-r 1 tan 2 §(fc + Z - j)u 



tan 2 JjV 



All further breather bound state intertwiners T are again uniquely given by the 
crossing relations using that the breather are self-conjugate 



r 



I 

k+lk 



r lk 
1 k+l 



1 lk+l 



j^kk-\-l 
1 I 



-pk+l 
1 kl 



-pk+l 
1 Ik ■ 



Again proposition ra and the general relation ftSA) imply that up to a sign 



Res Ski(9) 



Res S kk +i(9) 

=mr±(l-li/) 



which also may easily checked directly. 
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6.1.1 The bootstrap principle 

If there exist bound states in a quantum field theoretic model, the bootstrap 
principle means that all particles are to be considered on the same footing, in 
particular: 

1. The space of all kinds of particles V is closed under bound state fusion 
which means that there exist for each fusion angle a bound state inter- 
twiner \ u i2 ) such that 

lf 2 2) : Vi <g> V 2 -> V (12) , with ^1=^2 = V, V(i2) c V. 

2. If the fusion process a + (3 — > 7 exists then also the fusions + 7 — > a and 
7 + a — > /3 must exist. The corresponding fusion angles may be read off 
figure [?] where the euclidean momenta (p^Imp 1 ) are depicted and where 
u^p — 7r — etc. In particular the rapidities of the constituents are 
given in terms of the rapidity of the bound state by 

9 a = B 1 + iu^ 

e = ftf-iufo (40) 

and similar for the other fusion processes. The fusion angles satisfy the 
relation u 1 ^ + + u? a = 2tt or v? af3 + u^- + u^ a = n. 




Figure 7: The euclidian momenta (j>° jimp 1 ) of the fusion process a + (3 — > 7 
in the center of mass system of 7. 

3. The various bound state intertwiners are related by crossing as de- 
picted in figure || on page ||(| 

6.2 Bound state S-matrix 

We show that the 'bootstrap principle' provides a consistent scheme. In partic- 
ular we prove that the definition of the bound state intertwiners by the relation 
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34|) is consistent with crossing symmetry. First we define the bound state S- 
matrix (see p7j). Using the Yang-Baxter equation and the definition of the 
bound state intertwiners we have 



Res S12S13S23 — R es S23S13S12 
13S23 



r 21 r (12) r, 
1 (12) 1 12 °' 



— 593'5'l3r 





Therefore the following definition of the bound state S-matrix is natural. 



Definition 3 \2% 'The bound state S-matrix which describes the scattering of a 
bound state (12) with another particle 3 is given by 



5 , (i2)3(^(i2)3)r i 2 2) 



— T 12 S , i3(^i 3 )S , 23(023) 



1 ■ (12) 

> 12 =iu\ 2 



(41) 





where the rapidity of the bound state #(12] is defined by the relation of the 2- 
momenta p\ + P2 = P(i2) (see also eqs. HQ))- 

It was shown in [^7| that the bound state S-matrix defined by (Ej) fulfills 
unitarity, crossing and the Yang-Baxter equation. The equation (|4lf ) is also 
called a 'bootstrap equation' |3£j (also referred to as a 'pentagon equation'). It 
relates different two-particle S-matrices and therefore implies strong restrictions 
on the complete two-particle S-matrix which may be used to calculate the S- 
matrix for integrable models f^f . 

We show that the definition of the bound state intertwiners by means of 
the residue of the two-particle S-matrix (34) is consistent with the 'bootstrap 
principle' and crossing symmetry. 

Proposition 6 // the bound state intertwiners satisfy the crossing relations 



1 a/3 



X' 1 07'*-' 



(see figure and if the relation holds for the fusion process a + (3 
then it also holds for the fusion processes [5 + 7 — > a i.e. 



(42) 
-> 7 



i Res _ 



4f< 



_ p7'/3'p a 

— 1 a 1 /37 



and correspondingly for 7 + a — > [3. 

The proof of this proposition is delegated to appendix [d| 
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7 7 7 




a (3 a (3 a (3 



Figure 8: Crossing relations of the fusion intertwiners 



6.3 Bound state form factors 

If there are no bound states in the model there exist only the 'annihilation poles' 
according to property (Hi) on page |l^ which follow from the crossing formula. 
If there are also bound states there are additional poles jjj and we have also 
property (iv). Such additional poles have been discussed for simple cases in 
H Here we give the arguments for more general cases. Let us consider a 
model with a bound state (12) of two particles of 1 and 2 such that the attractive 
region is connected analytically (by a coupling constant) to a repulsive region, 
where the bound state decays. 

We start in the repulsive region and consider the two-point Wightman func- 
tion (JlJ> (0|O'O|0) = (0|C(x)O(0)|0). We use the symmetry given by the 
statistics of the particles to express the in-state matrix elements in terms of the 
form factor functions Oi,.. n (8) for ordered rapidities 



^ n\ J J (2^)«2w 1 . . 



dp x . . . dp n 

(27r)"2c i ;i . . . 2u> n 

x(0\O'\ Pl ,... ,p n ) mm ( Pn ,... , Pl \O\0) 
Mi f°° dd 2 f°° d6 n , nl 



n=0 



where P — pi + ■ ■ ■ + p n . For the last equality the symmetry (i) of the form 
factor functions has been used, in particular we have 

0' 12 ...n0 12 - n = 0' 2U „ n S i2 12 - n = 0' 1% ,J 21 21 - n . (43) 
In the repulsive case the S-matrices S\ 2 and S21 have poles at 812 = Q\ — 

(12) (12) 

62 = iu\ 2 and 6 21 = O2 — @i = iu\ 2 , respectively, in the 'unphysical region' 
U 12 < 0. Therefore also the left hand side has these poles. There are Q) such 

pairs of poles. By analytic continuation with respect to the coupling constant to 

(12) 

the attractive region where u 12 > these poles cross the integration contours 
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and we obtain additional contributions from residues 



(0 1 & o I 0) 



d9j_ 

47T 



47T 



0^ 
47T 



' ^12 


/•°° <%2) 


/•°° 1 


47T 


/-co 4^ 


7-oc 47T j 



-iPx 



The substitution (81,62) — > (6*12 = #i — #2, #(12)) nas been applied for the 
additional residue terms (if the masses of 1 and 2 are equal the bound state 
rapidity is 6(12) = \{9i +#2), for the general case see eqs. (4^)). Using (fl3|) and 
the residue formula (|34| ) we obtain 



(12) 47T 



12. ..rSS) - ~n 21...n(£) r (12) r i2 2) 



,(12) 47T 



1 r 21 r (12)^,21...n//i\ 

2 1 (12) 1 12 u 12-) 



] ■ (12) 

'l2=™i 2 



which means that both residues give the same contribution. Therefore the 
additional term may be written as 



E n! (2 

n=2 v 



(12) 



4tt 



dOn 
47T 



y p21 p(12)^-,l2...n 
'21. ..n 1 (12) 1 12 w 



1 ■ (12) 



dQ (12) 



E 



4tt ^(n-2)! 



d0 3 f 00 rf6>„ , (12) 
4^""./ 4^ U ( 12 )-" U 



iPx 



where we have introduced the bound state form factor 



O(12)3...n(0(12),0') = Ty| 21. ..n(£)r (12 ) 



] ■ (12) 

>12 = '» 12 



(44) 



with 6 = (63, . . . ,8 n ) and the rapidity 0(i 2 ) of the bound state is given by 
Pi + Vi = P(i2)- This identification is obviously unique up to a sign which may 
be absorbed into the bound state intertwiner. Using again the property (i) and 
the residue formula for the S-matrix ( |34| ) we obtain the property (iv) of form 
factors on page H 



Res Oi 

a - (12) 



23...n(£) = Res 2 i 3 ...„(£)5i2 

o (12) 

= -»0 21 3...»®r£ 2) rg? 2) 

012=P 12 

= 0(12)3...„(^12)^')V2r( 2 2) 

In order to show that the formula (^) for the bound state form factor is 
consistent with the bootstrap principle we prove the following proposition. 
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Proposition 7 The bound state form factor defined by the property (iv) or by 
(Hj satisfies properties (i) (ii) and (Hi) 

C(12)3...n(0(12): 03, ■ ■ ■ ) = ^3(12)...™ (03 , 0(12) > • • • )<%2)3 

= £ , 3...n(12)(03j • • • , 0(12) _ 27ri)a-0(i2) 
~ e (12 ^3-^ C (12)3C , 4...n(04,---) 

x (1 - S 3n . . . S 34 ) 

Proof. These relations follow directly from the corresponding relations for the 
form factor before taking the residue, by using in addition the Yang-Baxter 
equation, the fusion equation (fj~i| ) and the crossing relations of the bound state 
intcrtwiners. For (i) and (ii) the proofs are obvious. For (Hi) the proof is quite 
involved. It is delegated to appendix [5|. 



0) 
(ii) 

(m) 



7 Soliton Breather form factors 

In this section we apply the results of the previous section to the sine-Gordon 
model. We calculate breather form factors starting with the general formula 
( p"3| ) for the soliton form factors. 

The &fc-breather-(n — 2)-soliton form factor is obtained from Oi23... n (0) by 
means of the fusion procedure (iv) with the fusion angle given by u^P = u ^ — 
7r(l — kv) the bound state rapidity £ = &(i2) = |(0i + 02) an d = 03, • • • , n 

Res 123 ... n (6) = O(i 2) 3...„(£,0')v / 2r( 1 2 2) ( lM ( fc )) 

where the bound state intertwiner is given by eqs. (|3^) and (j39[). 

For 6\2 — + iu^ there will be pinchings of the integration contours in formula 
( [l3|) at the poles Zi — = Q^ — mlv = £ — |wr(l — kv+2lv) for / = 0, . . . , k and 
i = 1, . . . , to. Using the pinching rule of contour integrals and the symmetry 
with respect to the m z-integrations we obtain 

ReS 0123... n(0) 
#12— iu( k > 



After a lengthy calculation (see appendix [I]) we obtain for the case of the lowest 



k 

Res (— 2ixi)rnS Res / dz 2 ■ ■ ■ \ dz 

l2=i" (fc) ;=Q z 1= zW Jc L J CjL 



xh(e,z) P °(e,z)^i... n (e,z). 



38 



breather (k = 1, = 7r(l — the one-breather- (n — 2)-soliton form factor 



o 3 ...n(e,0') = n^K- fl< ) n 

2<i 2<i<j 
1 

X 



/ dz 2 --- dz m T[x(€- Zj,l) 

J Co J Co 



' Ce_ J Ce_ 

><nii^- z i) n '■(««)p o (e.£',« ( ' ) ,* / )*3...n(e',«') 

2<il<j Ki<j 

with z' = (z2, • • • j z m ). The soliton-breather form factor has been introduced as 

dt o cosh|^ 1 — coshi(l — 9 /(in)) 



i^hW = if s6 Wsinie exp / -2 

-cosfa-i/VBda-!/)) 



t cosh ^ 2 sinh £ 



sinhi(6»- f (l + v))ewhh(0 + ¥ (1 + v)) 



The normalization has been chosen such that F s b(oo) — 1. The function -E(i') 
was used in ^ ||] and is defined in appendix [f| Also we have introduced the 
short notations 

p(Cj0 = sinh|(e-(-l)'f(l + ^)) _ ( i); sinh|(e-f(l + (-l)M) 



sinh |£ sinh i£ 

,sinh±(£ + f (l + (-l)V)) 



sinh±(£-f (l + (-l)V)) 
The following identities have been used 

m-W2-^-^ (!) ) = F 8b ^-ei) P ^-e u i) 

^ - Zj )<f>(d 2 - Zj)S sb (Z ~ Zj ) T (z® - Zj ) = x (t - *j,t) 

for6» 1/2 = £±f (l-i/), = ^--f (l-(-l)'z/). The new p-function is obtained 
from the old one by 

P°(£,8', zW,z') = md{u) + - ±mW,0', Z W,£ ) 

where the constant d(v) is given by 



d{y) = 



y. x sin -,1:1' 



(see appendix |]). 
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Iterating the above procedure we obtain the r-brcather-s-soliton form fac- 
tor with 2r + s = n, the breather rapidities £ = (£1, • • ■ ,&•) an d the soliton 
rapidities 9 = (6i,...,6 s ) 

r s 

l<i<j'<r i= 1 j= 1 l<i<j'<s 



h=0 i r =0 l<i<j<r 
r s 



i sm 7Tf 
sinh ^ 



»=ij=i 

J i=l 7 =r+l 



s m 

i—lj—r+l r<i<j<m 



again with z^*' = £j — ^(1 — (— (i = 1, . . . , r). The two-breather form 
factor has been introduced as 



f dt 

F bb {9) = K bb {9) sini^exp / -2 



dt cosh(± - v)t 1 - coshi(l - 6»/(^tt)) 



o £ cosh it 2 sinh t 



cos T-irv/E{v) 



sinh i(6* — i7w) sinh ^(0 + iirv) 
The normalization has been chosen such that F bb (oo) = 1. The relation 

J^tfi - - 9 4 ) pfa - e s ,h)p(^ - 6i,h) X (Si - 4 h \h) 

„ / , j sin 7w 

= J J 66(€l2)(l + (il-i2)- 



sinh £i 

has been used for 3 / 4 = £2 ± -f-(l — f)- The new p-function is obtained from 
the old one by 

P(i,f,z^,z") = (™)r!fMp(( 1 |^^6-^ 1) ,.,«Vf l) ,.,/) • 
In particular for n = 2r = 2m we get the pure lowest breather form factor 

o(o = Il^(^) E-"E(-i)' i+ "" +,p 

i<j Zi=0 l r =0 



X 

l=*<j 
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and the pure breather p-function 

P% z!») = r! <f (v) p (fc + ±m« , & - ±m« . . . , z[ h) , . . . ) . 

8 Conclusion 

In a forthcoming paper p5| we investigate extensively the pure breather form 
factors of the sine-Gordon model. Some results have been published previously 
[ p4[ . Furthermore other integrable models of quantum field theory will be an- 
alyzed. In particular the S'C/(A^)-chiral-Gross-Neveu |l4|] and the (9(iV)-Gross- 
Neveu |l5| model is under investigation. In these models the particles possess 
anyonic statistics. The form factor program as considered in the present article 
may easily extended to anyonic fields and particles. Mainly the statistics factors 
in the form factor equations (i) — (v) have to be replaced by more general phase 
factors. This has been discussed in j|. [Io|. g, 
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Appendix 

A Proof of the lemma |2] 



Proof. We start the proof of lemma g by establishing some algebraic identities. 
The formula 



n ~ = 



l<i<j<m 



Xl 
X-2 



4 

■2 



1 X r 



„m — 1 



follows from the fact that both sides are polynomials in x m of degree m—1 with 
the same m—1 zeros and the same asymptotic behavior. Further taking the 
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anti-symmetric sum [•••]„ with respect to X\, . . . , x m one gets the relations 



Xm n x ^ 

l<i<j<m 



1 X\ 
1 X 2 

1 Xm 



m-2 k 



x m ^m 



= 8k,m-l Yi ( X l ~ 
l<.i<j<m 

for k = 0, 1, . . . , m — 1. Therefore for any set of constants one has 

/ rn — 1 \ 



.fe 



. fe=0 



n ^ ~ x ^> 

l<i<j<m 



l<i<j<m 



In particular with a?j = e 2zi and dj = e 26j we may finally write 



| j sinh z 

l<i<j<m 
'2m 

cosh \ ( 



'% **m 



) - ii sinn ^ 



i—l l<i<j<m 



because the difference on the right hand side is of the for m Y) T-n Ake 2kz ™ 



Using t{z) oc sinhzsinh(z/i/) the right hand side of formula (|l l| ) for p°(9,z) 
independent of z and n = 2m is proportional to 

/ dZi--- dz m Yi II ^ ~~ Z ^ I! Sillh ~ ( Zi ~~ Z ^ 



i=l j=l 



l<i<j<m 



J| cosh i(6>j - z m ) - J| sinh - z m ) j TJ e z J][ sinh; 

\i— 1 i—l / 1 l<i<j<m 

The z m -integral may be written as 

[ - ) dz m Y\_4>(0i - z m ) JJcosh|(^ - z m ) 



m — 1 



Y[ sinh ■ 



"I 



V(9,z) = 



where 0(0) = a(0)<f>(6) and *(0,z) = *(g, z)/ (n"=i HJLi a ( e > - «j) ) • The 
shift relations 

0(6* — i7w) cosh — (0 — i7ri/) = 0(0) i sinh -(0) 
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z u ... , z m + mv) = -$>(9, zx,...,z m ) 

have been used. They can easily be derived from (jl^), the definition of the 
Bethe Ansatz state (|l^) and the S-matrix ([!]) . The contour Cg — Cg_ + inv may 
be closed at ± infinity since the integrand behaves like e~ 2 ' Zm '/ l/ for zi — > ±00 
(see appendix |b|). There are no poles inside the closed contour, therefore the 
integral vanishes by Cauchy's theorem. ■ 



B Asymptotic formulae 

The asymptotic behavior of the soliton-soliton scattering amplitude is easily 
obtained from its integral representation 

a(9) = exp ^ g ffft 1 "?* riDhtm*)) = eVW"-* + o(l) 
Jo t sinh ^vt cosh 51 

for Re 6» -> ±00, (j Im6>| < f (1 + 1/ - |1 - This implies for the other 

amplitudes 

6(0) = e ±**(V"-i) + (i) , c (0) = ±2z sin(7r/^) e ±4(i/"-i) e T0/" (1 + ^ 

The asymptotic behavior of the 'minimal' two-soliton form factor function is 
given by 



r 

F(8) = cosh \ (i-n - 9) exp / 

Jo 



dt sinh|(l-i/)i 1 - cosht(l - 9 /(in)) 



Iq t sinh \vt cosh |t 2 sinh i 

= const e ±(1/v+1)(9 - i,r) (l + o(l)) 

for Re 6* — > ±00, (| Im0 — 7r| < -|(3 + ^ — |1 — v\)) with the constant 

dt ( sinh |(1 — v)t \ — v 



poo 

const = I exp ^ / 
Jo 



t V sinh ^z/t cosh ht sinh t i/t 



which satisfies 

const 4 = -\v (F'(0)) 2 = \vkm 

(for k see eq. ©). Similarly one has for Re(? — > ±00, (| lm8 — §| < §(2 + v 
|1 — i/| )) the asymptotic expansion 
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J ;e T|(i/H-l)(0-W2) 



Vn=0 
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for any integer n and v < 2/h. The constants A n are determined by the 
expansion 

dt ( sinhl(l-^)t _ 1-zA eit(e _ i7r/2)/7r " 
i \sinh i^i sinhi i^i / 

rn 

= exp(j2 a ™ eTnd + ( eTh9 ) + E ^ 2mdlv + o (e^™'") 

\n— 1 m— 1 

■n 

where a n = (— i)™" 1 sin \n(l — f)?V (7m sin \nvn} . Note that j4q = 1. 




C Expansion of the integral representation 

In order to compare the exact result for 4-particle form factors with the Feynman 
graph result the integral representation has to be expanded for small couplings. 
We calculate the integral 1^ = f c dz\ f c dz 2 I ± (zi, z 2 ) to prove formula (f28|). 
The integrands are 



i ± ( Z1 ,z 2 )= nn 



0(6*1 - Zi)c{V 2 - z 2 



v i=i j=i 



(l + b(6 1 - z 2 )b(6 2 - *!)) t{ Zi - z 2 ) (e ±Zl + e ±Z2 ) 



up to order 0(g 2 ). The functions I (zi, z 2 ) have poles at zi — 6\, #2, #3, 64 and 
9i — in, &i + in(u — 1) for i = 1, 2. By means of Cauchy's theorem we have 



dzil(zi, z 2 ) 



E 



dzi/(zi,z 2 ) 



where Co goes from —00 to 00 such that Im 9i — n< Im z < Im Oi . Then shifting 
by the integration contour Coby in one obtains 



dz i I(z 1 ,z 2 ) 



Co 



which implies finally 



E 

j= i - "i 



1) JCo+i-n 



dz i I(z 1 ,z 2 ) 



dzil(zi,z 2 ) = 



Co 



+I7r(l/-1) 



^E 



dzil(zi,z 2 ) + \ \ dz 2 (I(zi) + I(zi + in)) . 



C 
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The last term J c dzi (I(zi) + I(Z{ + in)) is of order 0(g 2 ). Therefore we may 
consider up to 0(g 2 ) 




The [— Res^ + Resg 1 -i n ] z . terms have O(l) and O(g) contributions the other 
ones are of order O(g) only. The 0(l)-terms cancel and after a long but straight 
forward calculation all 0(g)-terms give formula (|2g|). 



D Proof of proposition [6| 

In this appendix we prove that the 'bootstrap principle' provides a consistent 
scheme. In particular we prove proposition ^| which states that the definition of 
the bound state intertwiner is consistent with crossing symmetry. 
Proof. The idea of the proof may read off figure [|: We use the crossing property 
of the S-matrix (|J), the S-matrix bound state formula (|4l|), the residue formula 
( [34"| ) and the crossing relations of the bound state intertwiners (|35|). We make 




Figure 9: Proof of proposition 

the convention that rapidities as 9 a ,8 a ' etc. belong to particles with the same 
mass and other invariant quantum numbers. If S denotes a particle with the 
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same mass as that of j3 then we have (see fig. ^J) 

* Res a S}((0 M ) = i Res & C^' 7 'C 77 S^(wr - 9^) 

= i Res C^'T«?T^C^sfy(iir-6^) 

= -i Res c^'r^5f;«'(^)-5£ ,5 '(o)r^,,c 7 , 

S ofl=9 Q(3 

7' 7" a/3" a'/3^77 

/~<7'7'p/3'a'-p7 -p7'/3'-pa 

^ 1 7' i a'/3^77— 1 a 1 /3 7 

where the crossing relation of the S-matrix, the orthogonality relation of the 
bound state intertwiners (|37|), the residue formula ([mJ), the permutation prop- 
erty (||) of the S-matrix at = and the crossing relations of the intertwiners 
(HH) have been used. Obviously the rapidities satisfy the equivalence 

37 = iu % \ ^ / 0a/3 = *</3 



y « = tuls 



(note that # a /3 = 8 a — dp etc.). An analogous relation can be shown for the 
process 7 + a — > /3. The relation of the three fusion angles may be read off 
figure ■ 

E Proof of proposition [T] 

In this appendix we prove that the bound state form factors given by property 
(iv) on page || satisfy property (Hi) which is the recursion relation of n-particle 
form factors to (n — 2)-particle form factors. We consider (Hi) in the form 

Res 0(12)(34)5...n(#(12), 9(34), S.') = 2iC(i 2 )(34)05...ra(#') 

^(12)(34) — 

X (l - '5(34)n ■ ■ ■ <%4) 5 ) 

with ff = 9 5l . . . ,0 n . First we note that the bound state formula (iv) on page ffl 
for form factors may also be written in a form without taking a residue. This 
follows if we make the bound state pole at 9± 2 = iu explicit by using (i) such 
that 

C(12)3...n(#(12), = —7= Res 0123. ..„($!, 9 2 , ■ ■ ■ , 9 n )T}^ 2 

V 2 812=""' 



(12)' 



4= Res 213 ... n (9)S 12 (9 12 )T% 2) (-u) 

-i-^0 3 13...n(e) fl ^(u). (E.l) 
\/ L 9 12 =iu 



where also the definition (|34j) and the orthogonality relation (|37|) of the bound 
state intertwiners have been used. 



4G 



The form factor 01234...™ (0) has poles at 612 



^34 = IU, f 13 = «7T, 



724 



in. We extract the two first poles by applying formula (E.l) twice and introduce 
the function 



, r 21 r 43 
1 1 (12) L (34) 



G(0 12 ,0 34 ,8) — 0214:;:,...,, ! E I 1 ( L2) l. . \ti = P(i2)(34) = u (12) ~ °(34)) 

where we suppress the co- vector structure and the variables 64, . . . ,6 n . If for 
simplicity we assume that all masses are equal then 0(i 2 ) = \(Qi + #2) and 
0(34) = 5(6*3 + O4) are the center-of-mass rapidities. Extracting the poles at 
613 = in and #24 = in, the remainder is analytic and can be expanded as 

A + (6i2 - iu)B + (034 - iu)C + (0 - in)D 



G(012, 034, 0) — 



{013 - in) {024 - in) 
with certain constants A, B, G, D. We consider three limiting procedures 



(E.2) 



1. Let first 0i3 — > in, then $34 — > iu and finally — > in: 

This means that also $12 — > iu and 024 — > in. Due to (i) we have G(0i2, 034, 0) 

01324. ..n<S , 23<S f 2i5'43r^ L 2 )r^! 4 j and therefore with (Hi) for the particle pair 13 



d(0) 



Res G(0i2,0 3 4 

2i Ci3024...n (1 — S^n ■ ■ ■ S34S32) ^S^l^sF^^T 
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(12) x (34) 



2«Ci3024...nC43 ( — S^n ■ ■ ■ S35S21) ^12)^34) 



(2if 



024 - 

024 



-Ci3C24<T4 3 05...„(l — Sin ■ ■ ■ S45) ( — 5 3r . 



5rt \ y^21 t^43 
35<32lJ 1 (12) 1 (34) 



;Ci3C24<743<7i205... n (l _ S^n . . . S45) (— S3,, . . . S35) EJ^E^q 



~ — (2i) 2 «Ci 3 C 2 405...„ ( — S 3n . . . S 35 T^^ + r^| 4 - ) S'(34) n . . . £(34)5^ ^\i2) 

In deriving the third line unitarity (0) and crossing (JsJ) of the S-matrix has bee 
used. Also for (i, j) — (3,4) and (1,2) we have used that 

which follows from the unitarity of the S-matrix, the definition ([m|) and the 
orthogonality ( |37| ) of the bound state intertwiners. Further we have used again 
(Hi) for the particle pair 24 and the fact that (0i2 — iu)/ (024 — in) = — 1 holds 
for 0i3 = in and 034 = iu. The statistics factors 043012 cancel since 1 = 2 and 
3 = 4. Also the fusion rule S^SsiE^ = E^| 4 ^S , ( 3 4) i has been applied. On the 
other hand we obtain from eq. (E.2) using 024 — in = 2(0 — in) for 0i3 = in and 
012 — iu = —2(0 — in) for 0i3 = in and 034 = iu 

A — 2(0 — in)B + (0 — in)D 



d(0) 



.4 



2(0 -in) 



2(0 -in) 
-B 
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Since G\ (8) is non-singular for 8 — > in we conclude that A = and 

= -(2i) 2 zCi 3 C2405...ri . . . + r(34)S(34)n . . . £(34)5^ Tq 2 ) 

2. Let now #24 — * Mr, then #34 — ► zit and finally 8 — ► Z7r: 

Due to (i) we have G(8i2, 8 34 , 8) = 02413. ..riS^r^T^ and therefore similarly 
as above with (Hi) for the particle pair 24 



G 2 (8) = 



Res G(#i2,# 3 4, 



= 2i 



C24013...n (1 — 5*4™ • ■ ■ S45543541) "5*141^(12) T(34) 



2lC24Cl3...nO'21'S'2ir^ 2 - ) r^| 4 - ) 

(2i) 2 i- — Ci 3 C2405...n(l — Sa n . . . S 35)^} 12)^34) 

"l3 ^ i7T 

(2z) 2 «Ci3C24C5...ri(l ^ <?3n . . . "S^r^)! 1 ^ 



again we have used ( |E.3| ) and the fact that (#12 — iu)/(6\z — in) = 1 ho lds for 
6*24 = in and #34 = iu. On the other hand we obtain from eq. ( |E.2| ) using 
6*i3 — in — 2(8 — iir) for 024 = Mr and 6*12 — iu = 2(6* — i7r) for $13 = iir and 
6*34 = iu 



G 2 (8) 



A + 2{6-in)B- 



in)D 



2(8 -in) 



+ ... 



.4 



B+\D 



such that with A = 



G 2 (wr) = £+ |£> 

= (2 i ) 2 iC 1 3C 24 05...n(l-5; 



3ra ■ 



• (12) 1 (34) 



3. Let finally 12 = 6* 34 = and then 6* — > i7r: 

Applying formula (E.l) twice the form factor of two bound states and arbi- 
trary other particles is obtained as 0(i2)(34)5...n(#(i2) , 0(34), #') = — \G(iu, iu, 8). 
Therefore using 6*13 = 6*24 = 6 = 0(i2)(34) for 8\ 2 = 834 we obtain 



Res (12 ) (34)5...n (0(12), 0(34), £') = ~\ Res G(iu,iu,8) 

9(12)(34)=*f 9=47T 



\D 



\(G l (in)+G 2 (in)) 



2iC 13 C 2 A^\i2)^'(34)^ > 5...n(8')(l ~ 5*(34)n • ■ • 5(34)5) 
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which proves claim with the charge conjugation matrix for the bound states 



/~i -pl2 -p43 

(12)(34) - «^13^24i (12) 1 (34) 




(12) (34) 




(12) (34) 



F Residue formula 

In order to obtain the breather form factors from the soliton form factors we 
have to calculate the residue 

k 

Res 0i23 n (0) = Res (— 2m)m\^ Res / dz 2 ■ ■ ■ / dz m 

x h{e,z)p°{0,z)* 1 ... n (6,z) 
where the pinching rule of contour integrals has been applied. 
Lemma 8 For any analytic function P(z) the following identity holds 

Res Res Piz) MOx - ziW0 2 - z{) *i n (6,z) 

9 12 =iu( k > zi=z(0 

n m 

= c(k, l)(-l) l P(z^) J] a{6i - z<») J] S bkS (t - z j )lt 2 * 3 ... n (0 , ,z / ) 

i=3 j=2 

with & — 63, . . . , 9 n ,z? = Z2, ■ ■ ■ ,z m and 

*3...n(e',z') = n 3 ... n C 3 ... n (0', C 3 ...n(0', Z m ). 

The constant is 

-1/2 



c(k,l) = Rk-iRti 



\ Res S±(6) 



9=iw(l-kv) 



with 



Ri = Res 4>{z)a{z). 
Proof. Using the decomposition 

C\... n — C\ 2 A 3 n + -01203...™ 

and the action on the pseudo ground state f2 

Q12C12 = b{6i - z)c(6 2 - z) (s» s) 12 +c(0! - z)d(9 2 - z) (s® s) 12 
^121)12 = fii 2 6(0i - ^)6(6> 2 - z). 
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we obtain 

Res Res P(z) - z^O* - z x ) fii...„Ci...„(0,z) 

= P(«('))i2 fc _li2,(-1)' ((S ® S) 12 + (-l) k (s ® S) 12 ) 3 ...n^3. 

The relations 

Res = Res 0(6»)6(6») 

6) 12 =i7r(l-fejy) 9=t7r(l-(fe-Z)i/) 

= - Res <P(e)a(9) = -R k -t 

d—i7r(k—l)iy 

Res 0(#2 — z)d(# 2 — 2) = Res </>(— z)d(— z) = — Ri 

z=z(0 z=—i-Klv 



and 



c{mlv) _ j c(wr(l - (fc - 0^) 



a(iTtlv) 



b(in(l - (k - 



(_!)fc-J 

a(iTrlv) 



have been used. Further we use that the bound state intertwiner may be written 

as 



It 2 = (-l) fe i 



i Res 5±(0) 



1/2 



(s«)s + (-l) fe s®s) 12 . 



0=i7r(l-fci/) 

Applying further C-operators we use the fusion relation 

ri2 2) s lo (0 + ¥i2 ] )S2o{e - \e^) = s (12)0 (tf)rf 2 2) 









which implies 

r;$c7i... n (0, z)| eia=iuW - s 6fe5 (0 (12) - *)it 2 c 3 ...„(0', z) . 

Together with the eigenvalue equation 



n 3 ...„^3...„(£',2)=n^-*)°3.. 



j=3 



the claim follows. ■ 

Using the shift relation 



sin^(z + mv) 



(z)a(z) 
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the residues Ri can be calculated 



(— 1 V 

#i = #o- — i— r TT cot 5-? 7n/ 



sin i-7rZr; . 

2 3=1 

with 

a{z) —1 



i?o = Res (j)(z)a{z) — Res 



z=o^ v ' v ' 2=0 F(z)F(iir — z) F'(0)' 

In particular for the lowest breather k = 1 we obtain the constants independent 
ofZ 



c(l,0 = JVM 



-1/2 

i Res S±(0) 

z e=iir{l-kv) 



where is defined by eq. (19). Using these results the form factors for arbitrary 
numbers of lowest breathers and solitons are calculated explicitly in section ^. 
There we use the constant 

d{u) = ^Sf(mt(1 - v)) c(l, 0) = ^Vcos \vE{v) 



Sm 7T7V 



sm tjTT/' 



where the function E{y) is defined as 



50 di sinhz/i 1 tc/t 



-E(^) = CXp / = CXp — . 

Jo t 2cosh 2 ±i tt Jo smt 
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